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The locality problem for two tastes of staggered fermions * 
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We address the locality problem arising in simulations, which take the square root of the staggered fermion 
determinant as a Boltzmann weight to reduce the number of dynamical quark tastes from four to two. We study 
analytically and numerically the square root of the staggered fermion operator as a candidate to define a two taste 
theory from first principles. Although it has the correct weight, this operator is non-local in the continuum limit. 
Our work serves as a warning that fundamental properties of field theories might be violated when employing 
blindly the square root trick. The question, whether a local operator reproducing the square root of the staggered 
fermion determinant exists, is left open. 



1. INTRODUCTION 

The naive discretization of the Dirac operator 
on a Euclidean space-time lattice produces a the- 
ory with 16 Dirac fermion copies or tastes. The 
question whether this multiplicity can be reduced 
by a factor four has a positive answer. At finite 
lattice spacing the naive Dirac operator can be 
brought to a four-block diagonal form by spin di- 
agonalization. Each block is a local operator, the 
staggered Dirac operator, whose determinant is 
the fourth root of the determinant of the naive 
operator and describes four tastes in the contin- 
uum limit. The question which arises at this point 
is whether one can construct the "fourth root" (in 
the above sense) of the staggered operator. 

To this aim one should start from the taste ba- 
sis for staggered fermions. In the free theory four 
degenerate tastes can be identified in momentum 
space PP but are described by non-local fields in 
coordinate space and so this formulation is not 
suitable for introducing gauge fields. In coordi- 
nate space there is a basis of four taste-fields con- 
structed from paths which combine (in a gauge 
covariant way) staggered fields living in a 2^ hy- 
percube |2]- Even in the free theory there are 
taste-changing interactions, which vanish as 
(a is the lattice spacing) in the continuum limit. 
We would like to stress that the mixing terms 
of the hypercube fields are generated by a term 
with the structure like the Wilson term for Wil- 
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son fermions, which prevents the doubling of the 
four hypercube fields in the continuum limit. 

There is no obvious way of reducing the num- 
ber of staggered fermion tastes. For the Monte 
Carlo simulation of two degenerate (or one) dy- 
namical quarks one resorts to the so called square 
(or fourth) root trick P]. The Boltzmann weight 
is defined by taking the square (or fourth) root 
of the determinant of the staggered Dirac oper- 
ator, and this can be simulated efficiently. As 
concerns the observables, quark propagators are 
computed using the four-taste staggered operator 
and the sources project onto the desired valence 
taste components ^ . This construction seems to 
work in the classical continuum limit and in par- 
tially quenched chiral perturbation theory 5 . It 
is not based though on a first principle formula- 
tion of the fermionic theory, so far no acceptable 
operator describing two (or one) tastes has been 
found. Hence causality and the existence of a 
universal continuum limit of this theory cannot 
be established. 

One of the requirements for universality to hold 
is that the Dirac operator is local j6l7| . This 
means for its kernel in coordinate space G{x, y) 

\\G{x,y)\\ < ^7e-^ll-^'llB/^ ^>0, (1) 

where ||G(a;,j/)|| is the operator norm, ||a; — y\\E 
is the Euclidean norm and C, 7 are numbers in- 
dependent of the gauge field. Only local gauge 
paths are allowed in the definition of the Dirac 
operator. 
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In this work we present a locality test for a 
candidate operator to describe two tastes, namely 
the most naive choice, which is to take the square 
root of the staggered operator itself. This oper- 
ator has a determinant, which is obviously equal 
to the square root of the staggered determinant, 
but turns out to be non-local. The outcome 
might not be a surprise given the non-analyticity 
of the square root and the blindness concerning 
the taste structure. Nevertheless our investiga- 
tion can be considered as a benchmark study in 
the spirit of Ref. 6 for a still to be found oper- 
ator. Here wc only summarize our main results, 
which are published in detail in Ref. [H]. 

2. ANALYTIC STUDY OF A CANDI- 
DATE TWO-TASTE OPERATOR 

We denote by M the staggered Dirac operator. 
In numerical simulations the identity 

det(aM) = det{a^{M'' M)e) (2) 

is useful, where (M'l'M)e acts on fields defined on 
the even sites of the lattice only. We investigate 
the operator 



D = 



(3) 



where the square root is computed through a se- 
ries of Chebyshev polynomials. 
In the free theory 

AfU'f = rr? - A (2a) , (4) 

where m is the bare quark mass and A (2a) is 
the Laplace operator on a lattice with spacing 
2a. The kernel D{x) can be computed analyti- 
cally by Fourier transformation from momentum 
space. The result in the continuum limit is 
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(5) 



and means that the free operator D is non-local 
at the scale 

(free) 



1/to . 



(6) 



' loc 

It is very unlikely that the interacting theory will 
turn out to be local. At best one can hope that 
the localization range rjoc will be proportional to 
the inverse of a large hadron mass. This moti- 
vates our numerical study. 
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Figure 1. Decay of {f{r))/{f{rrei)), where r^ei is 
marked by the rightmost vertical dotted line. 

3. NUMERICAL RESULTS 

We work in the quenched approximation of 
QCD and generate gauge configurations at three 
values of the lattice spacing corresponding to 
/? = 6.0, 6.2, 6.5 with the Wilson gauge ac- 
tion. The bare quark mass m is tuned to keep 
the mass mc of the Goldstone pion ttq at the 
value roiriG = 1.30(3) (ro is the Sommer scale). 
We have two sets of volumes with lattice sizes 
LrriQ « 4 and Lttiq « 6. To study the locality 
of D we follow closely Ref. 6 . On each config- 
uration we define a source field S,c{x) to be 1 for 
one point x — y (we always set y = 0) and one 
color component c — 1, and zero otherwise. We 
compute the decay properties of ^(a;) = aDS^(x) 
in terms of 



/(r) = max{||i/'(a:;)|| I lla: - r} 



(7) 



Given a distance r from the source at y we find 
the maximum of the SU(3) color norm |j'0(a;)|| 
and denote it by f{r). For the distance r we take 
the "taxi-driver distance" 



y\\i = ^Tnm{\Xf, - y^l. 



\^t.-vA} (8) 



and its largest possible value is 2L. 

In Fig. n] we plot on a logarithmic scale the 
gauge average (/(r)) (normalized at the distance 
t-cci/to = 3.72) as a function of the distance r in 
units of tq. The results shown are for the volume 
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Figure 2. Effective localization range from {f{r)) 
computed at the two distances r marked by the 
vertical dotted lines in Fig. ^ 

Ltuq « 6 and the three lattice spacings. There 
is remarkable scaling as the continuum limit is 
approached, indeed suggesting a physical local- 
ization range. We assume the form 



(fir)) 



(9) 



and determine an inverse effective localization 
range rioc{r)~^^ using two consecutive values of 
(/(r)). The value of noc('') increases with r and 
has a "bump" at large distances r. This "bump" 
becomes higher as the physical lattice size L is 
increased and will eventually turn into a plateau 
on very large volumes, where there are no more 
finite volume effects. We understand the main 
origin of the latter 0. 

We adopt the strategy to compute the localiza- 
tion range at a fixed physical value of the distance 
r. For the choices r/ro = 1.86 and r/ro = 3.72 
(marked by the vertical dotted lines in Fig.^ we 
checked that there are no finite volume effects in 
noc('') computed in the volume Lmo ~ 6. The 
results for the different lattice spacings are shown 
in Fig. 121 and the continuum extrapolation yields 



riocir/ro = 1.86)tog 
riocir/ro = 3.72)tog 



0.455(12), 
1.06(5). 



(10) 
(11) 



^loc ™G 



2.8(6) 



(12) 



A local operator has r\oc{r) = 0(a) for all dis- 
tances r. The largest value of noc(^) in the vol- 



So we conclude that in the interacting case the 
scale of the inverse localization range is most 
likely given by the lightest pion mass. 

Our present work leaves the question open, 
whether a local operator D exists whose deter- 
minant is equal to the square root of the determi- 
nant of the staggered Dirac operator. We empha- 
size that present dynamical simulations, where 
the Boltzmann weight is taken to be the square 
root of the staggered determinant, do not use the 
square root of the staggered operator. If a lo- 
cal operator D is found then the configurations 
generated in these simulations are safe. There 
still remains however the unitarity problem. D 
will dictate the appropriate Green's functions for 
a two-taste theory and most likely these are not 
the ones built with the four-taste staggered oper- 
ator. 
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